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ABSTRACT 


The present work is aimed to explore combustion in a vertical diffuser with an emphasis on 
the buoyancy-induced nature of the flow. This combustion problem is similar to the 
combustion of pyrolysis gases in an annular vertical diffuser which is widely used in rural 
areas for cooking. The analysis of the problem is complicated owing to the coupling amongst 
natural convection flow, heat transfer and combustion. A simple finite reaction rate model for 
a single component fuel is presented for the combustion process. Then, the fuel is taken to be 
a mixture of combustible gases and an Arrhenius-type single step reaction is assumed between 
each of the combustible gases and oxygen. In each stage of modelling, temperature profiles 
are generated for various tube wall temperatures and fuel inflow rates. Results indicate that 
the maximum flame temperature is located close to the radial diffuser. The flame tends to 
move away from the wall with higher volatiles flow rate and higher wall temperatures also. It 
was also seen that maximum flame temperature increases with increase in tube wall 
temperature and power input, maximum flame temperature region spread in radial diffuser. 
Heat transfer efficiency increase with increase in power input. 


Keywords: Heat transfer efficiency, wood-volatiles, vertical diffuser, natural convection, 
mixture fraction formulation. 


1. INTRODUCTION 


Research and development of cook stoves and vertical combustors have been receiving its due 
share of attention at least in the developing world where a large chunk of the population still 
resides below a standard of living. This part of the population generally resorts to the use of 
wood, wood volatiles or any other bio mass fuel available easily for its cooking needs. The 
issue of stove or combustor is therefore quite sensitive and has direct implications on the 
socio-economic developments of the poorest of the poor human beings. A general layout of 
vertical diffuser combustor is shown in Figure 1 and 2. A combustion zone is sustained in the 
central passage by a buoyancy-induced upward flow of air and a radial inflow of combustible 
volatile gases from the vertical combustor. The hot flue gases transfer heat to the pan placed 
above the conical tube. Scientific investigation of such combustor is difficult because of the 
coupled phenomena of natural convection flow, heat transfer, and flaming combustion. In this 
work a certain radial flow rate of combustible volatiles has been prescribed from the inner 
surface of a vertical tube. In the present work, the combustion phenomenon of fuel gas in a 
buoyancy-induced airflow in a vertical combustor is modeled under the assumptions of a 
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finite reaction rate. The volatile flow rate and the inner wall temperature were parametrically 
varied over a restricted range for studying their effects on combustion. 
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2. MATHEMATICAL FORMULATION 


The flow, through the geometry shown in Figure 1, is modeled as axisymmetric. The 
computational domain is depicted in Figure 2. It consists of three regions, viz., straight tubular 
portion, and diffuser like conical portion known as axial diffuser and finally region between 
the stove and top disk is known as radial diffuser. The cylindrical polar coordinates r-z is used 
the present analysis. H is the total height of stove; s is the spacing between the stove and the 
pan. The flow is considered steady and turbulent and fluid is assumed to be incompressible. 
Since large variation of temperatures is encountered in stove, the Boussinesq approximation 
commonly used in natural convection flows is not employed in the present analysis. 
Consequently, all the fluid and thermodynamic properties including density are taken as 
functions of temperature. The equations of mass, momentum and energy describing the flow 
are written in non-dimensional form, represented by hat on the top of the variable. 


The present investigation mainly concentrates on modeling the actual combustion 
process along with fluid flow heat transfer. The fuel mixture is taken as combination of CHa, 
CoH6, CO, H2, CO2, and H20 which are typically part of volatile emerging from pyrolysis of 
wood. Composition of volatiles for CH4, C2H.e, CO, H2, CO2, and H20 is taken as 0.0489, 
0.10756, 0.4756, 0.0187, 0.249 and 0.1 respectively such that their combined calorific value 
becomes equal to that of typical wood sample. 


MASS CONSERVATION EQUATION: 


1 Oprv, OpV, 
oe ee (1) 
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RADIAL MOMENTUM EQUATION: 
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The symbols Gr, Pr and Da indicate the Grashoff number, the Prandtl number and the Darcy 
numbers respectively. Qsen refers to the source term representing uniform volumetric heat 
generation in the domain during combustion. 


SPECIES EQUATION: 


The combustion model consists of conservation of species equations. The fuel mixture 
contains species namely CH4, C2H«, CO2, H2, CO, HzO and Oz, No. These species are 
assumed to undergo a finite rate, single-step irreversible reaction kinetics. 

CH4+ 202——,CO2+ 2H20 

C2He + 3.5 O2 —, 2 CO2 + 3H20 

CO + 0.5 O2 ——, CO2 

Ho2+ 0.5 O2 ——, H20 
The equation conserving the species is then written as: 


lO7, a0) Ofaeasy 1 [le sof fe 
——{ prv Y. + V_Y. }= D D 0. 
sai (Nh) eb arose asl” (5) 


Le and D refer to the Lewis number and diffusion coefficient respectively while Y; and @j 
represent mass fraction and the rate of generation or depletion of the i™ species respectively. 
The rate of generation or depletion non-dimensional @; for each individual species mentioned 
above can be expressed in terms of Arrhenius equation. 
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3. BOUNDARY CONDITIONS 
3.1 AT THE FUEL INLET: 


The pressure at fuel inlet is taken equal to the stagnation pressure to the ambient pressure at 
the same elevation; the viscous losses due to acceleration of the fluid from surroundings to the 
stove inlet are neglected. Further, the fluid is considered entering the stove radially thus; axial 
velocity at the fuel inlet is taken to be zero. Then, continuity equation dictates the normal 
gradient of axial velocity is also zero. The temperature at the stove inlet can be assumed to the 
ambient temperature, 

aed < opV 


bos PV,» Vz=0, aie, 6=0 for O<?<1, 2=0 (6) 


3.2 AT THE AXIS OF SYMMETRY: 

Due to assumed axi-symmetry of the flow, computations may be carried out for only one axi- 
symmetric plane of the domain as shown in Fig. 1. Symmetry, then, dictates the radial 
velocity, normal gradients of the axial-velocity and the temperature respectively to be zero, 


. apv i i 
V. =0, .: =0, 0 for 7r=0, 0<Z<5+H (7) 
3.3 AT THE STOVES EXIT: 


The exhaust gases leaving the combustion chamber can be treated as a jet entering a still 
medium i.e. static pressure is equal to the local ambient pressure. Also, assuming that the flow 
leaves the exit radially, consequently, axial velocity is zero. Thus, boundary conditions at the 
exit can be written as; 


2 ane Q 
(Pq Py 82) . ép Vv. oT 
n= 7 V, = 0, ~ = 0, a FT 0 
P AG z or Or (8) 
3.4 AT SOLID WALLS: 


Ceramic liner is a solid wall boundary at all locations other than the primary and secondary 
air slots. At solid walls, the no slip boundary condition would apply, i.e., both radial and axial 
velocities are zero. For temperature, the convective boundary condition is imposed, 


: ; eT : -— 
V=0, V=0, rae ad for T= ax (2) 0<2<H (9) 


3.5 AT THE AIR INLETS: 

The air can be assumed to enter in the combustor axial direction and stagnation pressure is 
taken equal to the ambient pressure. The temperature at port can be taken as the ambient 
temperature. Air will flow due to free convection. 


1d 5 : 
p=- A.V, =0, =r =0, F=0 (10) 
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3.6 EQUATION FOR PROPERTY VARIATION: 
The fluid properties such as viscosity and thermal conductivity are assumed to be varying 
according to the Sutherland’s relation as follows 


3 ¢ 3/1 148 
» afl iS, | it , 19444 . 110.56 
eae p-#| | 2, fe —— 
6+S 6+S 
: ‘ (11) 


The variation of specific heat in terms of temperature is expressed by a cubic polynomial as: 
i a3 n2 4 
C, = 0.12917" + 0.34087” — 0.06917 + 0.9987 (12) 


3.7 EFFICIENCY OF HEAT TRANSFER: 


The efficiency of heat transfer is defined in terms of the net heat input to the gas as follows: 


Q an 
Nnt =o (13) 


Qgen—Qloss 
The heat transfer efficiency being an important result as far as performance of the combustor 
or stove is concerned. Since the variation in two quantities Qpan and Qioss depends on 
combustor geometry & volatiles flow rates. It has been seen in most cases that higher mass 
flow rates through the system leads to lower Qpan as well as lower Qioss due to lower 
temperatures. While a lower value of Qpan tends to decrease the efficiency of heat transfer 
lower Qioss tends to increase the efficiency. The net variation in nh.t. depends up on which 
quantity has more dominating influence. 


4. SOLUTION PROCEDURE 


A finite volume method was employed to discretize the governing equations. The SIMPLER 
algorithm, and standard «-e turbulence model were used. The discretized sets of algebraic 
equations were solved by Ansys—Fluent 13.0. To predict the effect near the wall clearly, 
uniform grids are selected in such a way that more number of grid points are near the wall. 


5. RESULTS AND DISCUSSION 


The variation in wall temperature and power input gives different combustion temperature. 
The temperature contours are generated of 800 K, 900K, and 1000K for 0.5 kW, 1.0 kW and 
1.2 kW for each case. The temperature profiles in the domain show that maximum 
temperature region moving away from the wall. This region also corresponds to the flame 
location in the combustion chamber. The maximum combustion temperatures for various 
cases are shown in Table 1. Maximum combustion temperature obtained is 1263 K. This 
observation also matches with maximum velocity region. 


Figure 4 shows variation of heat transfer efficiency for various input power. The simulations 
results are shown in Figure (a) to (i) for conditions given in Table 1. In general, these figures 
show that maximum combustion temperature increases with increase in wall temperature for 
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the same power input. It is clear that with increase in wall temperature flame region is 
increasing and flame is moving away from the wall and spread in radial diffuser. 


Similar effects are obtained with increase of power with constant wall temperature. It clearly 
indicates that maximum combustion temperature increases with increase in power and 
increase in wall temperature. 


The heat transfer efficiency being an important result as far as the performance of the 
combustor or stove is concerned. Since the variation in two quantities Qpan and Qioss depends 
on combustor geometry and volatiles flow rates. It has been seen in most of cases that higher 
mass flow rates through the system leads to lower Qpan as well as lower Qitoss due to lower 
temperatures. While a lower value of Qpan tends to decrease the efficiency of heat transfer 
lower Qjoss tends to increase the efficiency. The effect of variation in radial diffuser and stove 
power is depicted in Figure 4. The Qpan and Qioss are found to increase monotonically, so 
increase in power input combustion temperature increases, due to this fact heat transfer 
efficiency increases. 

The contours below show to temperature variation in vertical combustor. The contours 
arranged in the manner of case id according to Table 1.1. 


Table 1. Maximum combustion temperatures at glance for various simulations. 


Maximum Volatiles 
S.n. Case Id Wall Temperature Power Combustion Flow 
(KW) Temperature Velocity 
(K) (mm/s) 
1 PIT1 800 0.5 1190 5.20 
2 P2T1 800 1 1200 10.40 
3 P3T1 800 1.2 1206 12.48 
4 P1T2 900 0.5 1219 5.20 
5 P2T2 900 1 1230 10.40 
6 PSI 900 1.2 1232 12.48 
7 P1T3 1000 0.5 1249 5.20 
8 P23 1000 1 1262 10.40 
9 P3T3 1000 1.2 1263 12.48 
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Fig. 3. Iso-contours of temperatures in the computational domain for cases (a) P1T1, (b) P2T1, (c) 
P3T1, (d) PIT2, (e) P2T2, (f) P3T2, (g) P1T3, (h) P2T3, (1) P3T3 listed Tablel. 
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Fig.4. Variation of heat transfer efficiency with input power. 


6. CONCLUSION 


In the present work, the results of the combustion simulation are summarized as follows: 


(1) The temperature contours are predicted by the present combustion model. Results 
indicate movement of the flame away from the wall for higher volatiles flow rates and 
higher wall temperatures. 

(2) The region of maximum combustion temperature also corresponds the zone where the 
velocities are maximum. 
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(3) In general, maximum temperature in the flame region tends to increase with increase in 
wall temperature as well as power input. Maximum flame temperature region tends to 
move in radial diffuser. 

(4) Heat transfer efficiency monotonically increases with increase in power input. 
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